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SUMMARY 


A  method  is  developed  for  calculating  the  pressures  on  aerodynamic 
shapes  at  very  high  supersonic  speeds  in  dense  air  with  the  ratio  of 
specific  heats  equal  to  1.  The  method  is  applicable  to  any  body  1  of 
revolution  at  zero  angle  of  attack  and  to  any  two-dimensional  profile. 
The  results  of  the  present  paper  are  compared  with  previous  work  of 
von  Karman  and  Epstein  on  this  problem  and  the  differences  explained. 

Some  aerodynamic  characteristics  of  several  shapes  are  calculated, 
and  the  lift  and  drag  coefficients  are  shown  to  be  dependent  upon  the 
thickness  ratio,  thickness  distribution,  and  angle  of  attack. 


INTRODUCTION 


The  possibility  of  constructing  airplanes  and  missiles  capable  of 
traveling  at  extremely  high  supersonic  speeds  is  of  ever  increasing 
importance.  Reference  1,  for  instance,  shows  that  Mach  numbers  of  the 
order  of  30  or  bO  are  necessary  for  obtaining  very  long  ranges  with 
rockets.  The  largest  part  of  such  flight  paths  will  probably  be  at  very 
high  altitudes.  The  aerodynamics  of  flight  at  these  altitudes  has  been 
treated  by  Sanger  and  Tsien  in  references  2  and  3*  The  part  of  the 
flight  in  dense  air  (that  is,  in  a  fluid  which  may  be  considered  a  con¬ 
tinuum)  is  also  of  importance  since  the  surface  temperatures  and  air 
loads  may  reach  very  high  values.  The  problem  of  aerodynamic  heating 
has  been  investigated  in  reference  b .  It  is  the  purpose  of  the  present 
paper  to  investigate  the  air  loads  encountered. 

Busemann  has  done  much  of  the  original  work  in  developing  a  theory  for 
hypersonic  flows,  however  his  work  on  this  subject  is  not  generally  known 
in  this  country.  After  the  present  paper  had  been  completed.  Dr.  Busemann 
made  available  to  the  authors  a  copy  of  his  original  work  (reference  5)* 
The  method  developed  in  the  present  paper  is  in  substantial  agreement 
with  that  reference.  Von  Karman,  in  his  Volta  Congress  address  of  1935 
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(reference  6),  referred  to  Busemann's  theory  hut  in  applying  the  method 
imposed  physically  impossible  conditions  on  the  flow.  Another  theory 
has  been  presented  by  Epstein  (reference  l);  however,  since  he  neglected 
the  pressure  relief  due  to  centrifugal  force,  his  method  is  applicable 
only  to  wedges  and  cones. 


D 


F 

f 

L 

m 

M 

P 

P 

1 

r 

R 

t 

U 


*,y 


a 

3 

7 

P 


SYMBOLS 


drag  coefficient  based  on  unit  chord  or  maximum  cross  section 
lift  coefficient  based  on  unit  chord 


drag 


fineness  ratio 


Distance  from  nose  to  maximum  cross  section^ 
v  Maximum  thickness  ) 


function 

lift 


mass 

Mach  number 

pressure  coefficient 

static  pressure 

dynamic  pressure 

radius  of  body  of  revolution 

radius  of  curvature 

thickness  ratio 

free-Hstream  air  velocity 

coordinate  axes 

angle  of  attack 

angle  between  surface  and  free— stream  direction 

ratio  of  specific  heats 

density 
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0  shock-wave  angle 

0  angle  "between  surface  and.  x-axis 

Subscripts : 

o  free  stream;  due  to  form 

a,b  ahead  of  and  behind  shock  wave 

1,2  points  of  zero  pressure  coefficient 

c  due  to  centrifugal  force 

i  end  point  of  any  surface  element 

I  lower  surface 

L  due  to  lift 

max  maximum 

N  normal 

s  due  to  shock 

T  tangential;  total  due  to  pressure 

u  upper  surface 

METHOD  OF  ANALYSIS 


In  order  to  calculate  the  pressures  on  bodies  at  hypersonic  speeds 
it  is  necessary  to  determine  the  position  of  the  nose  shock  wave.  First 
consider  two  equations  from  the  theory  of  shock  waves.  The  relation 
between  the  density  ratio  across  a  shock  wave  and  the  Mach  number  normal 
to  the  shock  is 


7  +  1 


the  relation  between  density  ratio  and  the  shock-wave  angle 
deflection  angle  (3  is 

£b  _  tan  0 
PA  tan  (0  —  3) 


0  and 
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p 

If  Mxr  is  very  large  — —  can  be  neglected  and  the  following  very 

%2 

useful  approximation  for  high  Mach  number  flow  is  obtained: 


Pb 


7  +  1 
7-1 


For  7=1 


and  for  7=1.^ 


^->00 


—  — >6 

Since  9  is  not  generally  equal  to  90°  (tan  0  ^  °°) y  for  7  =  1 


e  =  p 

(la) 

7  =  1.4  and  slender  bodies 

0  =  1.2(3 

(lb) 

These  relations  may  he  applied  only  to  combinations  of  free— stream 
Mach  numbers  and  surface  angles  which  result  in  very  large  Mach  numbers 
normal  to  the  shock  (surface).  At  the  limit  M  =  oo  all  surface  angles 
may  be  considered. 

The  shock  equations  have  been  developed  on  the  assumption  of  an 
irreversible  adiabatic  process  (no  heat  transfer  except  that  occurring 
within  the  shock  wave  itself).  Epstein  (reference  7)  points  out  that 
the  compression  at  very  high  Mach  numbers  results  in  extremely  high- 
temperatures  which  in  turn  cause  large  heat  losses  due  to  radiation  and 
conduction.  These  heat  losses  limit  the  temperature  rise  to  a  value  much 
lower  than  that  determined  by  the  adiabatic  law.  A  better  approximation 
to  the  shock  equations  at  hypersonic  speeds  may  possibly  be  obtained  by 
use  of  the  assumption  that  7  =  1  in  the  preceding  equations.  Equa¬ 
tions  (l)  show  that  the  results  are  not  critically  affected  by  the  value 
chosen  for  7;  hence ,  for  simplicity ,  the  calculations  in  this  paper  will 
be  based  on  the  relation  7  =  1« 


Equations  (l)  indicate  that  at  very  high  supersonic  Mach  numbers 
the  shock  wave  follows  the  surface.  Since  the  flow  is  supersonic  the 
air  ahead  of  the  shock  wave  is  not  affected  by  the  body,  and  hence  the 
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field  of  disturbed  flow  around  the  body  is  limited  to  a  thin  layer  behind 
the  shock  which  has  been  called  the  hypersonic  boundary  layer.  The 
velocity  in  this  layer  may  be  determined  from  the  shock  equations. 

From  the  following  figure  it  is  apparent  that  the  mass  flow  (per 
unit  area)  through  the  shock  wave  is 

Pa^a^  “  Pb^bjj  (^ ) 


Since 


Pa 


Figure  1.—  Velocities  near  a  shock  wave. 
0  it  follows  that 


U- 


bN 


U, 


-*0 


aw 


Therefore  will  be  considered  negligible  relative  to  Ua^. 

Equating  momentum  before  and  after  the  shock  in  a  tangential 
direction  gives 


(paUaw)Uar  =  (p^N)UbT  (3) 

As  determined  from  figure  1,  equations  (l)^  (2),  and  (3)  the  tangential 
velocity  behind  the  shock  (along  the  surface)  is 

UbT  =  %  =  Ua  003  0  =  ua  cos  0 
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At  very  high  Mach  numbers  with  7=1  the  region  of  disturbed  flow 
is  confined  to  an  infinitesimal  layer  between  the  surface  and  the  shock 
wave.  Since  the  cross  section  of  this  disturbed  layer  is  very  small  and 
the  density  is  very  large,  the  acceleration  in  the  direction  of  flow  is 
negligible.  Hence  the  assumption  is  made  that  the  speed  of  a  given  mass 
of  the  fluid  remains  constant  along  the  surface  at  the  value  U0  cos  P. 
This  assumption  is  the  same  as  that  made  by  Busemann  and  von  Karman.  A 
sketch  of  the  velocities  over  a  surface  is  given  in  figure  2. 


cos  P  sin  Pj 


Uo  cos  P  cos  p^ 

^  Pi 


Airfoil  surface 


Figure  2.-  Velocities  over  a  surface. 

The  drag  of  the  body  is  found  from  the  rate  of  change  of  the 
momentum  of  the  air  in  the  stream  direction.  The  air  mass  per  unit  time 
entering  an  elemental  area  of  shock  for  a  body  of  revolution  at  zero 
angle  of  attack  is 

2itrp0U0  dr 

and  the  change  in  velocity  in  the  stream  direction  from  the  time  the 
elemental  mass  enters  the  shock  at  P  until  it  leaves  the  surface  being 
considered  at  ^  is 


Up,  —  Un  cos  P  cos  Pi 


The  drag  is  then 


2jcrp0U0(U0  -  UD  cos  P  cos  Pi)  dr 
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where  cos  p  is  a  known  function  of  r,  and  pp  corresponds  to  the 
surface  angle  at  the  point  rp  where  the  flow  diverges  from  the  surface. 

The  location  of  rp  will  be  determined  later  in  the  analysis.  Equation  (4) 
differs  from  that  given  by  von  Karman  inasmuch  as  he  assumes  that  the 
flow  leaves  the  surface  in  the  free— stream  direction  (Pp  =  0,  rp  = 

This  assumption  would  require  surface  pressures  less  than  absolute  zero 
and  is  therefore  not  physically  justified.  For  a  two-dimensional  body 
the  elemental  mass  is  p0U0  ^  ^ie  (^raS  ^or  one  surface  becomes 


D  = 


PoUo‘ 


(l  —  cos  p  cos  Pj )  dy 


(5) 


For  a  cone  P  remains  constant  at  the  same  value  as  Pp.  Then  the 
drag  of  a  cone  at  zero  angle  of  attack  simplifies  to 


max 


D  =  2 n  /  P0rU02(l  —  cos^Pp)  dr 

J0 


=  ^max^oV^in2^ 


(6) 


or 


C-n  =  2  sin^P- 


(7) 


Equation  (6)  is  identical  with  Epstein’s  result  (reference  7)*  Since 
Epstein  neglected  the  effect  of  the  pressure  drop  due  to  surface 
curvature,  his  method  would  be  expected  to  be  in  agreement  with  the 
present  work  only  for  straight— side  bodies. 

The  pressure  coefficient  behind  a  shock  wave  is  frequently  given  as 


2  sin  P  sin  6 
cos  (0  —  p) 


(8) 


Equations  (l)  and  (8)  combine  to  give  the  pressure  coefficient  behind  a 
shock  at  hypersonic  speeds 


P  =  2  sin2P 


(9) 
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A  comparison  of  equations  (7)  and  (9)  indicates  that,  for  7=1, 
two—  or  three-dimensional  todies  that  have  no  surface  curvature  in  the 
stream  direction  have  the  same  pressure  on  the  surface  as  behind  the 
leading  shock.  When  the  flow  is  curved  in  the  stream  direction,  the 
difference  in  pressure  from  the  shock  wave  to  the  surface  equals  the 
centrifugal  force  due  to  the  curvature  of  the  flow.  The  pressure  coeffi¬ 
cient  behind  the  shock  depends  on  the  surface  slope,  and  the .pressure 
relief  due  to  curvature  depends  on  the  local  air  mass,  velocity,  and 
radius  of  curvature.  The  maximum  curvature  that  the  flow  can  have  for 
a  given  shock  pressure,  mass,  and  velocity  occurs  when  the  pressure  on 
the  airfoil  side  of  the  flow  drops  to  zero.  If  a  pressure  still  lower 
is  needed  to  turn  the  flow  sufficiently  to  follow  the  airfoil  or  body 
surface,  it  is  physically  impossible  for  the  flow  to  follow  the  surface 
and  it  effectively  "separates."  For  the  very  high  dynamic  pressures 
encountered  at  hypersonic  speeds  the  pressure  coefficient  P  corresponding 
to  zero  pressure  can  be  considered  as  zero.  It  follows  that  the  limits 
of  integration  and  y-j_  of  equations  (4)  and  (5)  correspond  to  the 

points  where  the  surface  pressure  coefficients  go  to  zero.  Thus,  in 
order  to  determine  the  limits  of  integration  for  the  drag  integral  it  is 
necessary  to  know  the  local  pressure  coefficient. 

If  in  equation  (5)  the  subscripts  1  are  replaced  by  i  where 
the  subscript  i  corresponds  to  the  end  point  of  any  element,  the 
expression  will  give  the  drag  for  the  part  of  the  body  ahead  of  the 
point  (xqjyq) 


(11a) 
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and  for  xj_  ^  xq_ 


P±  =  0  (lib) 

Similarly,  for  a  three-dimensional  body  of  revolution  the  elemental 
drag  coefficient  based  on  the  maximum  cross  section  of  the  complete  body 
is 


Cj^  =  l6F2  J*  r(l  —  cos  0  cos  0^)  dr 
and  on  the  surface 

1 

?i  = 


For  xA  £  Xj_ 


o  2  sin  0<  d0^  rri 

2  sin  Pi  +  - - —  /  r  cos  0  dr 

ri 


(12a) 


and  for  Xj_  ^  x-^ 


P1  =  0  (12b) 

where  the  fineness  ratio  F  =  Length  from  floae  tcjmslam  c^os?  gec^on 

Maximum  thickness 

The  value  of  y^  =  7q  (or  r^  =  rq)  ie  determined  "by  equating  the 

pressure  coefficient  to  zero.  This  is  the  proper  value  to  be  used  in 
evaluating  the  total  drag  integral.  (See  equations  (4)  and  (5).)  The  flew 
separates  from  the  surface  at  the  point  (xq,yq)  and  is  bounded  by  a 
zero-pressure  streamline.  Any  body  -within  this  zero-pressure  region 
contributes  no  drag. 

In  equations  (lla)  and  (12a)  the  first  term  is  the  pressure  coeffi¬ 
cient  behind  the  shock  and  the  second  term  is  the  pressure  coefficient  due 
to  centrifugal  force. 
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A  good  approximation  to  the  surface  pressure  coefficient  for  slender 
profiles  may  he  obtained  in  a  simple  manner.  Let  y  =  f(x)  be  th^  shape 
of  the  body  or  airfoil  then 


and 


dy 

dx 


f’(x) 


f"(x) 


The  pressure  coefficient  behind,  the  shock  (equation  (9))  becomes 


ps  = 


2f ' 


1  +  f 


•2 


If  the  velocity  is  not  changed  appreciably  by  the  shock  (true  only  for 
slender  bodies  or  airfoils ),  the  pressure  drop  due  to  centrifugal 
force  is 


(Ap)c  = 


2  qy 

R 


or 


where  the  radius. of  curvature 


(1  +  f’-) 


3/2 


The  pressure  on  the  surface  is  the  shock  pressure  plus  the  change  due 
to  centrifugal  force: 


P  = 


Pa 


+  P, 


2f  *2  2ff 11 

1  +  f'2  (1  +  f '2)3/2 
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Since  f is  small  compared  with  1,  the  pressure  coefficient  on  a 
slender  airfoil  may  he  taken  as 


P  =  2(f'2  +  ff") 


(13) 


For  a  three-dimensional  "body  the  air  mass  over  each  surface  element 
is 


Po^2  pQr 
2nrU0  "  2 

By  similar  reasoning  the  pressure  coefficient  for  slender,  three- 
dimensional  todies  may  "be  approximated  as 

P  =  2f*2  +  ff"  (14) 

The  lift  of  a  hody  or  airfoil  results  from  the  downward  momentum 
imparted  to  the  air  stream.  The  same  considerations  used  in  determining 
the  drag  from  the  horizontal  momentum  apply  to  the  calculation  of  lift 
from  the  vertical  momentum.  The  lift  for  the  upper  surface  of  a  two- 
dimensional  airfoil  is  then  of  the  form 


L  =  -p0U0  sin  Pi 


*1 


cos  3  dy 


(15) 


where  y^_  is  the  point  where  the  surface  pressure  "becomes  negligible. 


APPLICATION  OF  THEORY 


Lift  and  drag  of  a  flat  plate.—  Equation  (9)  can  be  used  as  the 
basis  for  calculating  the  lift  of  straight— surface  airfoils.  Tne  lift 
coefficient  of  a  flat  plate  at  a  small  angle  of  attack  can  be  written  as 


Since  the  suction  pressures  are  negligible  in  comparison  with  the 
pressure  rises 


CL  =  2a2  -  0  =  2a2 


da 
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and 


Cp  =  C^a  =  2a3 

Lift  and  drag  of  wedge  airfoil.-  The  lift  of  a  single-^wedge  airfoil 
is  easily  obtained  by  the  same  method  used  in  the  preceding  section.  For 
angles  of  attack  a  less  than  the  semi-nose  angle 

CL  =  P2  “  pu 

CL  =  2  j{p  +  a)2  -  (p  -  a)2 


and 


CT  ~  8pa  =  4ta 


8p  =  kt 


(16) 


The  drag  coefficient  of  a  single-wedge  airfoil  at  an  angle  of 
attack  a  which  is  less  than  the  semi-nose  angle  is 

Cj^p  =  2  [(p  +  a)3  +  (P  -  a)3] 

=  4p3  +  12 Pa2 


This  can  be  broken  into  a  form  drag 


and  a  drag  due  to  lift 


=  12  Pa2 


(17) 


The  combination  of  equations  (l 6)  and  (17)  gives  the  following  expression 
for  drag  due  to  lift  for  this  airfoil: 
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For  angles  of  attack  larger  than  3,  the  lift  and  drag  of  the  wedge  and 
flat  plate  are  the  same  for  the  same  elopes  of  the  lower  surfaces. 

Lift,  drag,  and  pressure  on  a  two-dimensional  parabolic-arc  airfoil.— 
From  equation  (5)  the  drag  coefficient  for  a  symmetrical  two-dimensional 
airfoil  of  unit  chord  may  he  expressed  as 


C 


D 


(l  —  cos  3  cos  P}_)  d7 


For  the  parabolic  airfoil  defined  by  y  =  ±tx(l  —  where  t  is  the 
thickness  ratio,  the  drag  coefficient  becomes 


C 


D 


2tx1(2  -  xx) 


/T-Tt2 

1/1  +  t2d  -  Xl)2 


(18) 


The  value  of  x^  for  which  the  pressure  coefficient  is  zero  is  found 
to  he 


XX  =  1  -  i  /(I  +  t2)1/3  -  1 

Substitution  of  x^_  in  equation  (l8)  gives  the  following  equation  for 
drag  coefficient: 


CD  = 


2t  - 


(1  +  t2) 


pn1/3 


or  for  simplicity  of  numerical  calculations  it  may  he  expressed  in  a 
series  as 


+  •  • 
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Hie  pressure  coefficient;  on  the  surface  of  the  airfoil  at  any  point  (xi,yq) 
is  found  from  equations  (ll).  Then,  for  Xq  ^  xq. 


or,  approximately, 


P  =  t2  [3(1  -  Xi)2  -  1  +  1  +  6  (1  +  Xi)2  -  15(1  -  3Ci)^j  +  •  •  ' 

The  pressure  coefficient  found  by  the  approximate  method  (equation  (13)) 
is,  for  Xj_  ^  Xj, 


P  =  t2  [3(1  -  x^2  -  l] 

The  agreement  of  the  pressure  coefficients  found  by  these  two 
methods  is  very  good  as  can  be  seen  in  figure  3(a)  where  the  pressure 
coefficients  on  the  surface  as  well  as  the  pressure  behind  the  shock 
are  given  for  a  parabolic-arc  airfoil  of  thickness  ratio  equal  to  0.10. 

The  lift  coefficient  (based  on  a  unit  chord)  is  found  from 
equation  (I5)  to  be 

•  pJ2  r? 1 

CL  =  2  sin  p2  /  cos  0  dy  +  2  sin  Pq  /  cos  P  dy 

j  O  ^  O 
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where  (xi,yq)  (x2,yg)  are  '^■he  points  of  flow  separation  on  the 

upper  and  lower  surfaces,  respectively,  and  3  is  the  angle  between  the 
airfoil  surface  and  stream  direction.  Let  0  =  (3  —  a  for  the  upper 
surface  and  0  =  p  +  a  for  the  lower  surface  where  0  is  the  angle 
between  the  airfoil  surface  and  the  x-axis.  Then  for  small  angles 


sin  (0  —  a)  «»  0  —  a 

For  the  symmetrical  parabolic  airfoil  defined  by  y  =  tx^l  —  ^  the  lift 
coefficient  based  on  a  unit  chord  is  found  to  be 


=  [t^xg^  —  3txg^(t  +  a)  +  2xg(t  +  a)2J 

-Jt2xi3  -  3tX]_2(t  -  a)  +  2xq(t  —  a)2! 

The  pressure  coefficients  for  the  upper  and  lower  surfaces  are  found 
from  equations  (ll).  Setting  the  pressure  coefficients  found  from 
equations  (ll)  equal  to  zero  permits  the  determination  of  the  values 
of  Xj_  and  Xr> 


(19) 


Pu  =  3  t2x^2  —  6txj(t  —  a)  +  2(t  —  a)^  =  0 


P^  =  3"t2X22  —  6ti2(t  +  a)  +  2(t  +  a)2  =  0 


The  lift  is  then  found  by  substituting  these  values  of  x^_  and  Xq 
into  equation  (19).  For  a  <t 


(t  +  a)3  _  (t  -  a)^ 
t  t 


or  for  t  a  <  t 


1+1/3 


2/3 


CL  9 


(t  +  a)^ 
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or  for  a  '  t 


>♦.  (1  *<n 


CL  =  t2  -  3t(t  +  a)  +  2(t  +  a)' 


The  pressure  distribution  for  a  ten— percent-thick  parabolic  airfoil 
is  shown  in  figure  3(^)  for  a  =  0.05  radian.  If  a  is  larger  than  the 
semi-nose  angle,  the  effect  of  the  upper  surface  on  the  total  lift  is 
zero  since  suction  pressures  are  neglected. 

Drag  and  pressure  on  a  parabolic  body  of  revolution.—  The  drag  of  a 
parabolic  body  of  revolution  at  zero  angle  of  attack  may  be  determined  by 
the  same  method  used  for  the  two-dimensional  bodies.  Let  the  surface  be 
described  by  the  equation 


-  K1  - 1) 


where  F  is  the  fineness  ratio  of  the  body 
coefficient  is  then 


The  expression  for  drag 


ri 


CD  =  16F2  /  (1  —  cos  0  cos  cLr 


which  reduces  to 


CD  =  2x^(2  -  Xlr  + 


2  8F 


If2  +  1 


+  2  - 


(1  -  x-|  )2  (F2  +  1) 


F2 


-  2 


F2 


'1  + 


(l  ~  Xj)' 
F2 


J 


where  x^  is  again  the  point  of  zero  pressure  coefficient. 
The  expression  for  Cp  expanded  in  a  series  is 

1  !1  +  3(1  -  x;l)2  ~  9(1  -  x-l)’4  +  5(1  -  xx)6 


=  3^  L’ 


24f 


3  +  4(1  -  Xi)2  +  18(1  -  xx)^  -  60(1  -  xj_f  +  35(1  -  xx) 
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The  surface  pressure  coefficient  found  from  equation  (12)  is,  for  xq  ^.xq. 


or  as  a  series,  for  xq  <  Xq, 


12Xq(2  -  Xq ) 


|l  -  6(1  -  Xq)' 


(20) 


By  us©  of  the  approximation  given  by  equation  (l4)  the  pressure  coefficient 
may  also  he  expressed,  for  x^  <  x^_,  as 


P  =  ^2  (5xi2-10xi  +  V  (21) 


The  pressure  distribution  over  a  parabolic  body  of  revolution  of  fine¬ 
ness  ratio  equal  to  10  as  determined  by  equation  (20 )  is  plotted  in 
figure  4.  The  value  of  the  pressure  distribution  as  determined  by  the 
approximate  method  (equation  (21 ))  is  also  plotted  in  figure  4  together 
with  the  pressure  behind  the  shock. 
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CONCLUDING  REMARKS 


A  method  for  calculating  the  aerodynamic  characteristics  of  two— 
and  three-dimensional  shapes  at  hypersonic  speeds  in  dense  air  has  been 
developed.  The  results  of  the  present  paper  are  in  substantial  agreement  ^ 
with  some  earlier  work  of  Busemann  but  differ  from  the  results  of  von  Karmen 
and  Epstein.  Von  Karman,  in  his  application  of  Busemann* s  method,  has 
assumed  that  the  flow  leaves  the  surface  parallel  to  the  free  stream. 

This  assumption  would  require  surface  pressures  less  than  absolute  zero. 
Epstein  neglected  the  pressure  relief  due  to  centrifugal  force  and  there¬ 
fore  his  results  apply  only  to  bodies  with  zero  curvature. 

The  method  of  the  present  paper  has  been  applied  to  several  simple 
shapes  and  compared  with  other  calculations. 


Langley  Memorial  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  February  5,  19^8 
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(a)  a  =  0  radian;  Cp  =  0.00066. 

Figure  3.-  Pressure  distribution  over  parabolic  airfoil.  Thickness  ratio 


20 


NACA  TN  No.  l6l3 


a 


Pressure  distribution  over  parabolic  body  of  revolution.  Fineness  ratio,  10;  a  =  0  radian. 


